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Abstract. The coronal slow shock has been predicted to exist embedded in 
large coronal holes at 4-10 solar radii. We use a three-fluid model to study 
the jumps in minor ion properties across the coronal slow shock. We formulate 
the jump conditions in the de Hof f mann-Teller frame of reference. The Rankine- 
Hugoniot solution determines the MHD flow and the magnetic field across the 
shocks. For each minor ion species, the fluid equations for the conservation 
of mass, momentum and energy can be solved to determine the velocity and the 
temperature of the ions across the shock. We also obtain a similarity solution 
for heavy ions. The results show that on the downstream side of the coronal 
slow shock the ion temperatures are nearly proportional to the ion masses for 
He, 0, Si, and Fe in agreement with observed ion temperatures in the inner 
solar wind. This indicates that the possibly existing coronal slow shock can 
be responsible for the observed heating of minor ions in the solar wind. 
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1. Introduction . 


Alpha particles and other minor ions have been extensively observed since 
the first measurement of the solar wind. The dynamical behavior of alpha 
particles and other minor ions in the solar wind is not yet well understood. 
The helium to hydrogen abundance ratio e - n a /n^ is highly variable. The long 
term average value of £ is in the range of 0.03 to 0.06. The average value of 
the speed ratio, U a to Up, is slightly greater than 1.00. The speed ratio, U a 
to Up, near 0.3 AU is greater than that near 1 AU. The most probable value of 
T a /Tp is between 3 and 4. There is a strong tendency in the solar wind for the 
ion temperatures to be roughly proportional to the masses [Bochsler, 1989; 
Bochsler, et al., 1985; Hernandez, et al., 1987; Marsch, et al., 1982; 
Neugebauer, 1981; Neugebauer and Snyder, 1966; Ogilvie, et. al., 1968; 

Ogilvie, et al., 1989]. 

In this paper we study the jumps in flow properties for alpha particles and 
other minor ions across slow shocks. A few slow shocks have been identified in 
the solar wind [Richter, 1987], but the changes in properties of minor ions 
across them have not been published. Observations of changes of alpha 
particles velocity, density, and temperature have been published, both for 
fast interplanetary shocks and for the earth's bow shock [Borodkova et al., 
1989; Fuselier at al., 1988; Neugebauer, 1970; Ogilvie et al., 1982; Zastenker 
et al., 1986]. Both Ogilvie et al. [1982] and Bovodkova et al. [1989] reported 
that reacting to the presence of the cross-shock potential difference the 
alpha particles are decelerated less than the protons and a temperature-mass 
proportionality is introduced by the action of a fast shock. It has also been 
observed that magnetic forces contribute appreciably to the slowing of ions in 
the de Hof f mann-Teller frame of reference [Thomsen et al., 1987]. We expected 
that slow shocks and fast shocks may have similar effects on the dynamical # 
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behavior of solar wind ions. This study concludes that the possibly existing 
coronal slow shock can be responsible for the heating of minor ions and 
perhaps the variation in flow velocities for different kinds of ions in the 
solar wind. 

The coronal slow shock has been predicted to exist embedded in large 
coronal holes at 4-10 solar radii [Whang, 1982]. The model considers the 
heliomagnetic polar regions of the sun as the sources of solar wind streams, 
interplanetary magnetic fields of one polarity emanating from north polar cap 
and fields of opposite polarity from south polar cap are separated by a 
neutral sheet. In coronal space, the high-speed solar wind streams emanating 
from the polar coronal holes are sub-Alfvenic flows of low-/9 plasma. Studies 
for the expansion of the solar wind from coronal holes in rapidly diverging 
stream tubes indicate that supersonic speeds should be attained at points low 
in the corona [Kopp and Holzer, 1976; Munro and Jackson, 1977; Whang, 1983; 
Whang and Chien, 1978], Streams originating from the edge of the polar open- 
field regions flow around the curved boundary of the helmet- shaped, closed- 
field region. At the edge of the neutral sheet, the flow direction of each 
stream changes suddenly, becoming parallel to the neutral sheet. Due to 
dynamical interaction in north-south direction between neighboring stream 
tubes, an oblique MHD slow shock can develop near the neutral point. The shock 
is upstream inclined, and extends polewards to form a standing coronal slow 
shock surrounding the sun. 

We use a three- fluid model to study the jumps in the flow properties of 
minor ions across slow shocks. We formulate the jump conditions in the de 
Hoffmann-Teller frame of reference. The Rankine-Hugoniot solution determines 
the MHD flow and the magnetic field across the shocks. The main body of the 
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paper presents a three fluid model to calculate the flow of protons and alpha 

particles across slow shocks under the assumption (a) n Q is a small fraction 

of , but not negligible in comparison with, n^ and (b) the flow properties of 

electrons across slow shocks is specified by a polytropic relation. The model 

1 6 6 + 28 8 + 

is also used to calculate the flow of other minor ions 0 , Si and 

^Fe^ + ) across slow shocks under the limits of n^ « V 

The fluid equations for the conservation of mass, momentum and energy for 
each kind of ion can be solved to determine the velocity and the temperature 
of the ions across the shock. If the protons and each kind of minor ions (k - 
He, 0, Si, and Fe) have the same temperatures and the same velocities on the 
upstream side of a slow shock, then the model shows that under reasonable 
coronal conditions, on the downstream side of the coronal slow shock the 
calculated ion temperatures are nearly proportional to the ion masses and the 
average ion velocities are slightly greater than 1. Assuming that these 
conditions are preserved during the expansion of the solar wind between the 
corona and 0.3 AU, the consistent results between this study and the 
observations of minor ions in the inner solar wind suggest that these 
observations support the theoretical prediction for the existence of the 
coronal slow shock. We also obtain a similarity solution showing that on the 
downstream side of a standing slow shock the temperatures are directly 
proportional to masses for heavy ions. 
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2. Three-Fluid Model 


2A. MHD Shocks 

The classical Rankine-Hugoniot Relation is used to study shocks in a MHD 
fluid consisting of electrons, protons and alpha particles. We choose a 
cartesian coordinates system: the xy-plane is the plane of coplanarity, the x 
axis is normal to the shock surface pointing the direction of the mass flow, 
and the tangential component of the magnetic field points in the positive 
direction of the y-axis. We formulate the flow conditions in the de Hoffmann- 
Teller frame of reference. Outside the shock layer, the MHD fluid velocity is 
align with the magnetic field B, and the electric field E and the non- 
coplanarity component of the magnetic field B z are zeros. Inside the shock 
layer E - Ee x , an< * E, and the fluid properties for each kind of particles 
(k - e, p and a) are functions of x only. This formulation includes a 
noncoplanarity component of the magnetic field inside the shock layer 
[Goodrich and Scudder, 1984; Gosling, et al. p 1988; Jones and Ellison, 1987; 
Thomsen, et al M 1987]. The equations for the conservation of mass, momentum 
and energy are used to study the change in flow velocity and temperature for 
minor ions across the shock. 

The jumps in the magnetic field B and the fluid properties p, U, and T 
outside the shock layer can be determined by the flow conditions upstream of 
the shock using the Rankine -Hugoniot relation: 

[pU x ] - 0 (1) 

pU x (U] + [P + B 2 /8ff]e x - B x [B]/4* - 0 (2) 

[U 2 /2 + c p T ] - 0 (3) 

and 
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[B x ] - 0 


(4) 


Here the pair of square brackets denote the jump of a physical quantity across 
the shock 


[Q] - Q 2 - Qi < 5 > 

the subscripts 1 and 2 respectively denote the flow conditions upstream and 
downstream of the shock. 

Whang [1987] has a simple direct method to calculate the solution of the 
Rankine-Hugoniot relation to determine the jumps in the MHD flow and the 
magnetic field across oblique MHD shocks. We use this method to calculate the 

ratios B 2 /B 1# p 2 /p ^^xl’ P 2/ P l and the shock an S le 9 2* < The shock an S le 9 
is the angle between the shock normal and the magnetic field). They are 
functions of three dimensionless upstream parameters: the shock Alfven number 
A — u x l/( a l cos 0]_)» the shock angle 9 and the plasma ft value. (Here a is 
the Alfven speed, and 0 the ratio of the thermal pressure to the magnetic 
pressure A shock is a slow shock if the shock Alfven number A less 

than 1 and is a fast shock if A is greater than 1. The magnetic field and the 
shock angle decrease across a slow shock, they increase across a fast shock. 

Shock solutions have been systematically studied in a three dimensional 
parametric space [Edmiston and Kennel, 1986; Whang, 1988]. For given /3 and 9 ^ , 
shock solutions exist in the domain of A > A min . Here the minimum shock Alfven 
number 

A min 2 “ U + 5/3/6 - (d + 50/6) 2 - 10 / 9cos 2 fl 1 /3} 1/2 )/2cos 2 ^ 1 (6) 

for slow shocks, and 

A min 2 “ U + 5£/6 + { (1 + 50/6) 2 - lQ0cos 2 e l /2) 1/2 )/2cos 2 8 1 (7) 


8 



I 


for fast shocks. Evolutionary shocks do not exist in the domain o f A < A rain- 
There are no jumps in MHD flow properties and magnetic field when A - A m ^ n> At 
any given value of 0 , a whole range of changes in from 0 to 1, may take 

place in the solution domain. The range of change in thermodynamic properties 
strongly depends on the 0 value. Slow shocks covering a wide range of jumps in 
thermodynamic properties (such as p and p) exist in the low 0 region (0 < 

0.1). On the other hand, only slow shocks with weak jumps in thermodynamic 
properties exist in the high 0 region (0 > 1). Therefore, in the coronal space 
where the 0 value of the solar wind is of the order of 0.1, all physical 
properties may jump across a slow shock over a wide range of magnitudes. Near 
1 AU the plasma 0 value is in the order of 1, the jumps in thermodynamic 
properties across a slow shock must be small but the change in the magnetic 
field is not necessarily small. 

In place of the shock Alfven number, shock Mach number is sometimes used as 
an independent variable to calculate the shock relationships. The shock Mach 
number is defined as the ration of U n to the magnetoacoustic speed, C g for 
slow shocks and C f for fast shocks. In the M,/9,0 parameter space, slow shock 
solutions exist in the domain of 

1 < M 2 < (1.2 + 0 + ((1.2 + 0) 2 - 4 . 80 cos 2 0) 1/2 )/2/3 (8) 

2B. Three Species of the MHD Fluid 

We denote the three species of the MHD fluid, electrons, protons and alpha 
particles, respectively by subscripts e, p and a. We consider that the number 
density of alpha particles n a is a small fraction of, but not negligible in 
comparison with the proton number density np. Inside the shock layer the 
charge density is in the order of (h/5)^n e e, where h is the Debye shielding 
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distance and 8 measures the shock thickness. Since h « 8, we may assume that 


to the zeroth order of magnitude 


n e “ "p + 2n a (9) 

In the de Hoffmann-Teller frame of reference, the electric field inside the 
shock layer is generated entirely by the charge separation which is a first 
order quantity. Now we can express the density and the flow velocity of the 
fluid consisting of three kinds of charged particles as 


P " “p 1 ^ 1 + 4 0 


( 10 ) 


U + 4€U 
P « 

1 + 4f 


(ID 


Since this study is concerned with different velocities and temperatures 
for alpha particle an protons, we should introduce the diffusion velocities 


and 


W 


U 


U 



( 12 ) 


They are respectively the velocities of alpha particles and protons in a frame 
of reference moving at the fluid velocity U. is expected to be very small 
compared with W . Depending on the frame of references, there are two pressure 
tensors for alpha particles: (P *)j^ the ii component of the pressure tensor 

in the frame of reference moving at the fluid velocity U, and (P a )^ in the 
frame of reference moving at the flow velocity of alpha particle U a , 


^ii “ ^a^ii + m a n a W ai 


( 13 ) 


A similar relationship exists between (P and (Pp)^- The fluid pressure 
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calculated from the Rankine-Hugoniot equations is the average of pressure 
tensor in the fluid frame of reference 


P 



2 

i 


{ (Oii + <P« >M + CP« >M ) 


p 'ii 


’ ii 


( 14 ) 


On the other hand, the temperatures of alpha particles and protons are 
respectively defined in the frame of reference moving at flow velocities of 
alpha particles and protons. Namely 


and 


nkT - 
a a 


— S < P a>ii 
3 i 


n_kT_ - 
P P 


— 2 ( p p)ii 
3 i K 


(15) 


From these relations, we obtain that the fluid pressure 

P - ^(npTp + n a T Q + n g T e ) + + n> a n a W a 2 ) (16) 

In the next two subsections, we will develop a three- fluid model to 
calculate the flow of alpha particles across slow shocks. The equations for 
the conservation of mass, momentum and energy are used to calculate the change 
in n , U a and T Q . Once we have the solutions for the alpha particles and for 
the MHD fluid, then we can calculate n^, T^, and U^. 

2C. Lorentz Forces 

Inside the shock layer, the electromagnetic field exerts a Lorentz force on 
each kind of charge particles (k - e, p, a) 

- Z k n k eE + i J k x B (17) 

where Z g - -1, Zp - 1, Z Q - 2, and e is the elementary charge. This force and 
its work done on the moving plasma affect the momentum flux and the energy 
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flux of the charged particles across the shock layer. Note that the parallel 
component of the electric current density makes no contribution to the Lorentz 
forces. We may express the Lorentz force as the sura of two parts 

H - h? * V < 18 > 

The first part represents the electric Lorentz force due to the electric 

field and the electric drift current density and second part is the 

magnetic Lorentz force due to perpendicular electric current caused by the 
sudden change in the magnetic field configuration inside the shock layer. 

We can express the Lorentz force due to the electric field and the electric 
drift current density as 

L k E - Z k n k e {E + (E x B) x B/B 2 } 

or C 

- Zk n k e ^ cos ^ e l (19) 

where e-^ » B/B is the unit vector along the field direction. Summing over all 
kinds of charged particles, 

2 L k E - 0 (20) 

k 

This means that the electric Lorentz force produces no net effect on jumps of 
fluid properties for the MHD fluid, the mixture of all kinds of charged 
particles, as shown in the Rankine-Hugoniot equations. 

Making use of Ampere's law, the sum over all kinds of charged particles of 
the Lorentz forces caused by the magnetic field configuration inside the shock 
layer must be equal to (V x B) x B/47T 

2 L k M - (V x B) x B/4 n (21) 
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IntegraCion of this equation over the shock layer produces the accumulated 
effects of the Lorentz forces which appears in the Rankine-Hugoniot equation 
(2). Equation (21) provides an important constraint the magnetic Lorentz force 

. However, it is difficult to formulate the exact expression for each . 

u 

Our scheme is using the guiding-center theory to estimate the ratios L e : Lp 
: L , then can obtain a reasonably representation for each of them from (21). 

The perpendicular electric current densities driven by the magnetic field 
configuration can be calculated from the motion of the guiding centers which 
are the instantaneous centers of the particle orbits. (The parallel component 
of the electric current density makes no contribution to the Lorentz forces.) 
Using the gyration theory of charged particles we can describes the effect of 
the changing magnetic field configuration on the motion of the guiding center. 
For each kind of charged particle there are three electric current densities 
related to the magnetic field configuration: the gradient drift current, the 
curvature drift current and the magnetization current. From these drift 
current densities, we can obtain the portion of driven by the magnetic 
field configuration as 


j ik 


< p iik - p ik> e i x * 


( 22 ) 


where k - e-^.Ve^ is the curvature vector, || and J_ refer to directions relative 
to the local magnetic field, ic points the instantaneous center of the circular 
arc and the radius of curvature equals to 1/k. This equation shows that the 
perpendicular electric current density J ^ (k - p, a, e) inside the shock 
layer is directly proportional to the thermal pressure anisotropy (P^ * Pj^). 
Since the electrons are nearly thermally isotropic, our model assumes that 
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(23) 


L. M - 0 

Based on this formulation, we can obtain an estimate for the 

h" ■■ l « m - < P IIP - P 1P> : < p n° - V 

From Equations (21), (23) and (24) we obtain 

L ^ 1 - r k (V X B) x B/4 w 


following ratios 

(24) 

(25) 


for k - p and a where 


r 


r 


p 


f hp ~ p ip 

(p hp ■ p ip ) + (p ii« * P i a) 

p l|q ' p l<* 

(p hp - p ip ) + (p n« - V 


(26) 


(27) 


2D. Jump Conditions for Alpha Particles 

The equations for conservation of mass, momentum and energy for alpha 
particles can be integrated to give 


t n a U aX l - 0 

t n a m a U ax 2 + n a kT cJ “ F a 
f n a m ct^ax^ay ^ 

and 

[U a 2 /2 + 5kT a /2m a ] - - Z a e[*]/m a 

where 


(28) 


(29) 


(30) 


(31) 
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( 32 ) 


f x d B 

F - { Z ft n eE — - - I\ — (— ) ) dx 

a J a * R 2 ° dx 8ir 


- | < Z o", 


B B B dB 

«e 2LZ + r Q ji -Z> d* 

a ^2 a 4tt dx 


(33) 


and <£ is the electrostatic potential. The term on the right hand side of 
Equation (31) represents an energy sink. The alpha particles do work as they 
move against the electrostatic field inside the shock layer. F a and G a are 
respectively the integration of L ax and L a y across the shock layer. They 
represent the accumulated effects of the Lorentz forces over the shock layer 
on the momentum flux of alpha particles. The first parts of F a and G a shows 
that alpha particles are decelerated by the electric Lorentz forces. The 
second parts of F Q and G a represents the effects of the magnetic Lorentz 
forces on alpha particles. Inside a slow shock, the magnetic Lorentz forces 
accelerate alpha particles in the shock normal direction and decelerate alpha 
particles in the tangential direction. Slow shocks occur in low f} plasma such 
as in the solar coronal or in the geomagnetic tail environment. Equations (32) 
and (33) show that except for parallel shocks the magnetic Lorentz forces 
normally play an more important role than the electric Lorentz forces in 
affecting the dynamical behavior of alpha particles crossing slow shocks. The 
integrations for the z component of the momentum equations are not shown here 
because they only produce first order corrections to the flow speeds. 


Equations (28)-(31) are used to calculate the ratios n aX 2/ n a xl’ U ax2 //U axl f 
U ay2 /U a yi and ^ a 2/^ a l across the shock layer. We can first use (31) to 
calculate the jump in tangential velocity across the shock 


IV 1 


G 


a 


m a n al U axl 


(34) 
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Then we solve Equations (28), (29) and (31) for n QX 2/ n a xl’ ^ax2/ U axl an< * 
T a2 /T Q ^. The system of equations can be reduced to a quadratic equation for 

^ax2^axl 


< U ax2/ U axl)' ' ^axZ^axl) + *0 “ 0 


(35) 


where 


5 kT al 

A 1 “ — { 1 + — + 


m a U axl 2 m a n al U axl 2 


(36) 


A 0 - — 


i tV ] 


5kT 


al Z a e W 


(37) 


4U axl 4ra a U axl' 


2m U 


a axl 


The equation has two roots associated with the + signs of the square root 
term. The solution with - sign represents the shock solution. This step 
produces the most important dynamical effect of shock wave on alpha particles 
part of the kinetic energy is converted to increase the thermal energy. Once 
we have U ax2 /U axl , we can calculate n a2 /n Ql , T a2 /T al from 


and 


n a2/ n al “ U axl/ U ax2 


(38) 


T a2 “ T al * (m Q /5k)[U Q 2 ] - (2Z a e/5k)[*] 


(39) 
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3. Results for Alpha Particles 


3A. Computational Model 

Section 2 provides a theoretical formulation for the jumps in flow 
properties of protons and alpha particles across slow shocks. The same 
formulation can be used to calculate the jump of other minor ions across slow 
shocks (in section 4) and to calculate the jump of ions across fast shocks (in 
a separate paper). In order to carry out numerical computations, we need to 
make three further assumptions with regard to (a) the flow conditions upstream 
of the shock, (b) the electron flow across the shock, and (c) the thermal 
anisotropy of protons and alpha particles. Note that once we have better 
understanding about electron flow and thermal anisotropy of ions, the 
theoretical formulation obtained in this paper can be upgraded to carry out 
new calculations. 

We assume that on the upstream side of the shock all kinds of particles 
have the same flow velocity and the same fluid temperature 


T al 

" T pl “ T el 

(40) 

U al 

- U pl “ U el 

(41) 


We assume that inside the shock layer the thermodynamic properties of 
electrons may be described by a polytropic relation 

P e " P el< n e/ n el > 7 

(42) 

or 

T e - Tel *"./".!) 1 ' 1 

where 7 is the polytropic index ranging between 5/3 and 3 [Feldman, 1985; 
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Schwartz, et al., 1987; Schwartz, 1988; Scudder, 1987; Thomsen, et al., 1987]. 
Computations carried out in this paper use 7 - 5/3. Although the electrons 
carry a significant parallel electric current, the parallel electric current 
makes no contribution to the Lorentz forces. According to the gyration theory 
of charged particles, the Lorentz force induced by the magnetic field 
configuration on electrons is negligible because the electrons are nearly 
isotropic [Equation (23)]. In the equation of motion for electrons, the 
inertia term is negligible because the electron kinetic energy is much less 
than the electron thermal energy. The x-component of the equation of motion 
for electrons can be reduced to 

dP e 2 

- - n^eE cos z 9 (43) 

dx e 

Under the polytropic law one can calculate the electrostatic potential rise 
across a shock from (43) 

e [^ ] - 7 k<cos' 2 0 >[T e ]/( 7 - 1). (44) 

Here the pair of angle brackets denote the average value of a physical 
quantity inside the shock layer. 

We assume that the distribution functions of the solar wind alpha particles 
and protons have similar asymmetries with respect to the magnetic field line. 
The alpha particle thermal anisotropy is generally slightly, but perhaps not 
significantly lower than that for protons near 1 AU [Bosqued et al., 1977; 
Hundhausen et al., 1970; Marsch et al. 1982; Ogilvie et al., 1980]. The ratios 
p j| p/ p j^p an< * [a Can f unct ^ ons of x. If we assume that ^j|p/^J_p “ 

P|| a /Pj^ a , then we can write (26) and (27) as 
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(45) 


r p“ 


P__ + P„ 
P a 


r a " 


(46) 


P„ + P„ 
P a 


Now, we can express the integrals for F a and G a in the following form: 


"a <»y 2 l 

- - 2. < > [f«l - <r > 


(47) 


n. 


8tt 


and 


n cr B y 

G o - - Z o < ~ > [ p e ] + <r a> ~ 


(48) 


n e B x 


4tt 


Since we do not know the detailed structure of the flow and field inside the 
shock layer, we have to use an approximate method to estimate those average 
values of physical properties which appear in equations (44), (47) and (48). 
We use the average of the values on the two sides of the shock layer to 
represent the average, namely <Q> - (Q^ + Q 2 V 2 . If Q 2 is an unknown variable 
of the equation system, an iteration scheme is used to calculate the average 
values . 

We calculate the flow conditions downstream of slow shocks, with a 
particular interest in the two ratios T a 2 /Tp 2 > an< * ^a2^p2* ^ or a S^ ven 
combination of the plasma f} value, the shock angle 0^, and a dimensionless 
shock speed, the shock Mach number or the shock Alfven number A - U x ^/(a^ cos 
0^), we first calculate the ratios ^x2 //U xl > U y2 // ^yl f ^2^1 and 

02 using the simple direct method of Whang [1987] for the MHD Rankine -Hugoniot 
relation. Then we calculate the flow conditions of alpha particles for the 
ratios n a2 /n al> U ax2 /U axl , U ay2 /U Qyl and T a2 /T al across the shock layer 
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following the method developed in section 2D for a given the initial helium to 
hydrogen abundance ratio " Hal^Hpl* ^ird ste P is to calculate the flow 

conditions of protons using formulas developed in section 2C. The solutions 
can be organized as functions of four dimensionless parameters: y9, plus 

A or M. 

3B. Alpha Particles 

The shock relations formulated in the de Hof fmann-Teller frame of reference 
have a singularity at 9 ^ — 90°. The iteration scheme also becomes difficult to 
converge for nearly perpendicular shocks. We carry out numerical solutions for 
0 ° < < 60°. For a given combination of 0 and 6^, we can construct constant 

contour plots for ^ a 2/^p2 and ^a2 / ^ U p2 on M »^l pi ane as shown in Figures 1 and 
2. The thick lines represent the constant contours for T^/T 2 anc * tllin li nes 
for U Q 2 /Up 2 * The two panels in Figure 1 show that the variation in the helium 
to hydrogen abundance ratio, - 0.05 and 0.10, produces very insignificant 
changes for T a 2/Tp2 and U Q 2 /^p 2 * The domain of solution for slow shocks is 
bounded on the right by a curve representing M - M max . Slow shocks are 
evolutionary for 1 < M < M max . Four contour plots with varying yS values are 
shown in Figure 2. The plasma fi value has a stronger effect on the solutions 
of T a 2/Tp 2 and U a2 //U p2- The varia tion in P substantially changes the upper 
limit of the slow Mach number. As changes from 0.01 to 0.20, M max decreases 

from -11.0 to =s2.5. At very small M m ax var *- es as 7(l-2/0). 

The flow parameters upstream of the coronal slow shock may be estimated to 
be in the range of y9 « 0.02, M > 3.0 and 15° < 9^ < 45° [Whang, 1982, Figure 
1]. From Figure 2, we find that immediately downstream of the coronal slow 
shock the temperature ratio T a /T p should be in the range between 3 to 4 and 
the average speed ratio U a /U p should be slightly greater than 1. The ion 
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temperatures approximately proportional to the ion masses can be expected to 
preserve between the corona and 0.3 AU. 
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4. Other Minor Ions 


4A. Computations 

The three-fluid model can be used to calculate the conditions of other 
minor ions across a slow shock. In this case the MHD fluid consists of 
electrons e, protons p and a species of minor ions i - ^®Si® + , ^Fe^ + 

or others. We can assume that n^ « rip. Under this assumption, we can use n^ - 
n 0 , Up - U, and T ^ - <P^/Pp> in the computation. Figures 3-5 shows the contour 
plots for the ratios to the fluid velocity U 2 and to fluid temperature 
T 2 respectively for and ^Fe^ + at two values. Once again 

these figures show that the temperature ratio proportional to ion mass is the 
most important dynamical effect of shock wave on minor ions. The coronal slow 
shock convert a part of the kinetic energy to increase the thermal energy. A 
careful examination of Figures 4 and 5 reveals a striking result that the 
constant contours for and T-^/ 111 ! are a l raost: identical in these plots. 

4B. Similarity Solutions for Heavy Ions 

For heavy minor ions, because is much greater than 1 Equations (32) and 
(34) show that the magnetic Lorentz force becomes less important as compared 
with the electric Lorentz force. Because m^ is much greater than nip, in the 
upstream side of the slow shock the thermal energy is much less than the 
kinetic energy. If we neglect these small order terms, then we can arrange the 
jump conditions for heavy ions in the following form: 


n i2 

— U ix2 “ U xl 
n il 


(49) 


n i2 2 kT i2 

— < U ix2 + 

il 


z i f n i B ' 
— eE 
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> - u x r + 
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m. 


dx 


(50) 
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dx 


(51) 


n 12 Z i ( n i B x B y 

— U ix2 U iy2 - U xl U yl + “ ~ eE 

n il m i J n il ® 

and 

U 12 2 5 kT 12 u l 2 Z 1 

+ e[<f>] (52) 

2 2 ra^ 2 

Here we have a system of four equations for four unknowns n i2^ n il r ^ix2* U iy2» 
and kT^/®^- In th® limit of n^ « protons and electrons are responsible 
for the setting up of the electric field inside the shock layer and the 
variation of the magnetic field across the shock layer. If we consider two 
kinds of heavy ions having the same ratios Z^/m^ across the same shock wave, 
the two flows are governed by the same system of four equations. All 
corresponding terms on the right hand sides of the two systems are identical. 
The two equation systems have the same solutions. The two flows are described 
by a set of similarity solutions. This similarity relationship explains that 
no noticeable differences exist between the computational results in Figure 4 
for Si and in Figure 5 for Fe. The similarity solution for kT^/m^ means that 
across the coronal slow shock the temperatures are proportional to the masses 
for heavy ions. 
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5. New Support for the Coronal Slow Shock 


We use the fluid equations for the conservation of mass, momentum and 
energy to describe the ion flow across a postulated coronal slow shock 
predicted to occur in the low (3 coronal space between 4 and 10 solar radii. 
From figures 1 - 5 we can see that on the downstream side of the coronal slow 
shock, for He, 0, Si, and Fe the ion temperatures are approximately 
proportional to the ion masses and the average ion velocities are slightly 
greater than the proton velocity. The result is consistent with the ion 
temperatures and velocities observed in the inner solar wind. We infer this 
result as a new support for the possible existence of coronal slow shocks. 

If the coronal slow shock indeed exists, our understanding of the solar 
wind process will be significantly revised. There exists one illustrative 
example which calculates the solar wind flow associated with the possibly 
existing coronal slow shocks [Whang, 1986]. The example has calculate a strong 
latitudinal variation of the solar wind, the terminal speed at the pole is 
greater than that at the equator by a factor of 2 and the plasma density at 
the pole is less than that at the equator by a factor of three. The low- 
density and high-speed solar wind at the pole region should be directly 
observable from Ulysses spacecraft. 

The support provided by this study of minor ions across slow shocks should 
stimulate modelers to consider the coronal slow shock in their solar wind 
models. A systematic study of coronal slow shocks should produce more 
quantitative solutions from which we can estimate the range of variation for 
important shock parameters such as the position and the shock strength and 
identify the range of coronal conditions under which the coronal slow shocks 
may exist. Because the predicted coronal slow shock is within the reach of the 
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future Solar Probe spacecraft, continued theoretical studies of the coronal 
slow shock will become very helpful in the planning of the Solar Probe mission 
and its experiments. 
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Figure Captions 


Fig. 1 The thick lines represent the constant contours for ^ a 2/^p2 anc * c ^ n 
lines for ^ a 2/^p2’ two P ane ^ s show that the variation in produces very 

insignificant changes for ^ a 2/^p2 anc * ^ J a2 / ^p2‘ 

Fig. 2. Four contour plots for slow shocks with e ^ - 0.05 and - 0.01, 0.02, 
0.10 and 0.20 show that the plasma j 3 value has a strong effect on the 
solutions of T a 2 /Tp 2 and 

Fig. 3. The three-fluid model is used to calculate the conditions of other 
minor ions across a slow shock. Two contour plots show the ratios V^/^2 and 
T i2 / T 2 for i - 16 o 6+ with 0 - 0.01 and 0.1. 

Fig. 4. Two contour plots for i - with /? - 0.01 and 0.1. The ion 

temperature ratio proportional to ion mass is the most important dynamical 
effect of shock wave on minor ions. 

Fig. 5. Two contour plots for i - ^Fe^ + with - 0.01 and 0.1. We also 
obtain a similarity solution for heavy ions to explain the reason that no 
noticeable differences exist between the computational results in Figure 4 for 
Si and in Figure 5 for Fe. 
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Fig. 1 The thick lines represent the constant contours for ^ a 2/^p2 an< ^ c ^in 
lines for U a2 / U n2- The two P ane l s show that the variation in ^ produces very 
insignificant changes for T Q 2 /Tp 2 anc * ^a2^p2' 


30 




0.01 e= 0.05 T./T, U„/U p /3 = 0.10 c=0.05 

on on u n vn fi. G 7.0 8.0 9.0 10.0 11.0 12.0 1.0 2.0 



31 


Fig. 2. Four contour plots for slow shocks with 6-^ - 0.05 and = 0 . 
0.10 and 0.20 show that the plasma value has a strong effect on the 
solut ions of T 0 /T. r, and H )/U >. 
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Fig. 3. The three-fluid model is used to calculate the conditions of other 
minor ions across a slow shock. Two contour plots show the ratios U^/Uo and 
T i2 /T 2 for 1 " 16 ° 6+ with £ - 0.01 and 0.1. 
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Fig. 4. Two contour plots for i - 28 Si 8+ with 0 - 0.01 and 0.1. The ion 
temperature ratio proportional to ion mass is the most important dynamical 
effect of shock wave on minor ions. 
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0=0.01 i= 56 Fe 16+ T,/Tp U , /U p 
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Fig. 5. Two contour plots for i - 6 Fe* with p - 
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